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Frustrated Electron Liquids in the Hubbard Model* 
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The ground state of the Hubbard model is studied within the constrained Hilbert space where 
no order parameter exists. The self-energy of electrons is decomposed into the single-site and 
multisite self-energies. The calculation of the single-site self-energy is mapped to a problem of 
self-consistently determining and solving the Anderson model. When an electron reservoir is 
explicitly considered, it is proved that the single-site self-energy is that of a normal Fermi liquid 
even if the multisite self-energy is anomalous. Thus, the ground state is a normal Fermi liquid 
in the supreme single-site approximation (S 3 A). In the strong-coupling regime, the Fermi liquid 
is stabilized by the Kondo effect in the S 3 A and is further stabilized by the Fock-type term of 
the superexchange interaction or the resonating-valence-bond (RVB) mechanism beyond the 
S 3 A. The stabilized Fermi liquid is frustrated as much as an RVB spin liquid in the Heisenberg 
model. It is a relevant unperturbed state that can be used to study a normal or anomalous 
Fermi liquid and an ordered state in the whole Hilbert space by Kondo lattice theory. Even if 
higher-order multisite terms than the Fock-type term are considered, the ground state cannot 
be a Mott insulator. It can be merely a gapless semiconductor even if the multisite self-energy 
is so anomalous that it is divergent at the chemical potential. A Mott insulator is only possible 
as a high temperature phase. 

KEYWORDS: Mott insulator, metal-insulator transition, Hubbard model, single-site approximation, Kondo 
effect, Kondo lattice, Fermi liquid, RVB, frustration, third law of thermodynamics 



1. Introduction 

It is an important but longstanding issue to elucidate 
the nature of the Mott metal- insulator transition. ^ Al- 
though many studies have been performed, the issue is 
still contentious. The third law of thermodynamics ap- 
pears to be broken in a Mott insulator, as discussed be- 
low in this section. One of the most contentious issues is 
whether or not the ground state can be really a highly 
degenerate insulator or the third law can be really bro- 
ken in a physically relevant model, or the Hubbard model 
with a finite on-site repulsion U, which is one of the sim- 
plest effective Hamiltonians to study this issue. 

According to Hubbard's theory, 2, 3 ) when the on-site 
U is as large as or larger than the bandwidth W, a band 
splits into two subbands; the subbands and a gap be- 
tween them are called the upper Hubbard band (UHB), 
the lower Hubbard band (LHB), and the Hubbard gap, 
respectively. Hubbard's theory can be pictured by a sim- 
ple argument. First, consider the case of U/W — +00. 
The number of electrons per unit cell is denoted by n. 
When n = 1, all the unit cells are singly occupied and 
there is no empty or double occupancy. The half-filled 
ground state is the highly degenerate state of having an 
electron of arbitrary spin at each unit cell; the entropy 
is fee In 2 per unit cell. This ground state is a prototypic 
Mott insulator. Assume that an electron is removed from 
the Mott insulator or a hole is introduced to it. There is 
no reason why the hole is localized in a periodic system, 
or the Hubbard model; the bandwidth of the hole mo- 
tion is 0(W), i.e., the bandwidth of LHB is 0(W). When 
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U/W = +00, the ground state is an insulator for n = 1 
and is a metal for 1. Next, reduce U. LHB is centered 
at the original band center e a and UHB is centered at 
e a + U. Their bandwidths are modified by the virtual ex- 
change processes allowing empty and double occupancies 
but are still 0(W). When U is sufficiently large, the Hub- 
bard gap cannot close. Thus, one may argue that when 
n = 1 and U is sufficiently large the ground state is an in- 
sulator. However, this argument is too naive to conclude 
that, when U is finite, the ground state can be an insu- 
lator. According to Gutzwiller's theory 4-6 ) together with 
the Fermi-liquid theory, 7 ' 8 ) a narrow quasi-particle band 
forms near the chemical potential; the band is called the 
Gutzwiller band in this paper. The ground state can only 
be an insulator if the Gutzwiller band disappears. Ac- 
cording to Brinkman and Rice's theory, 9 ) which is in 
the Gutzwiller approximation, the Gutzwiller band dis- 
appears when n — 1 and U > Ubr, where Ubr — W . 
In the Gutzwiller approximation, the ground state is an 
insulator when n = 1 and U > Ubr and is a metal when 
n 7^ 1 or U < Ubr- The divergence of the specific-heat 
coefficient as U — > Ubr implies that the third law ap- 
pears to be broken in the insulator. 9 ) 

One may speculate based on both Hubbard's and 
Gutzwiller's theories that the density of states has a 
three-peak structure, with the Gutzwiller band between 
UHB and LHB. The approximations used in Hubbard's 
and Gutzwiller's theories are within the single-site ap- 
proximation (SSA); they are called the Hubbard and 
Gutzwiller approximations, respectively. According to 
another SSA theory based on a slave-boson or auxiliary- 
particle Hubbard model, 10 ) the Gutzwiller band forms at 
the top of LHB when n < 1, which implies that it forms 
at the bottom of UHB when n > 1. The spatial dimcn- 



2 J. Phys. Soc. Jpn. 



Full Paper 



F. J. Ohkawa and T. Toyama 



sionality is denoted by d. The SSA that considers all the 
single-site terms is rigorous in the limit of d — > +00, n ~ 14 ) 
but within the constrained Hilbert space where no or- 
der parameter exists; 15 - 1 this SSA is called the supreme 
SSA (S 3 A). The S 3 A theory is mapped to a problem of 
self-consistently determining and solving the Anderson 
model, 16 ~ 18 ) which is an effective Hamiltonian to study 
the Kondo effect. The three-peak structure corresponds 
to that in the Anderson model, with the Kondo peak 
between two subpeaks. The Kondo effect is crucial for 
electron correlations in the Hubbard model. The S 3 A is 
also formulated as the dynamical mean-field theory 19-22 ) 
(DMFT) and the dynamical coherent potential approxi- 
mation 23 ) (DCPA). 

It is interesting to develop \ jd expansion or Kondo lat- 
tice theory, in which leading order terms in 1/d or single- 
site terms are non-perturbatively considered within the 
S 3 A, DMFT, or DCPA and then higher-order terms in 
1/d or multisite terms are perturbatively considered be- 
yond it. In an early Kondo lattice theory, 16 ~ 18 ) which 
was proposed about a year earlier than the first paper 
of DMFT, 19 ) it was assumed that the ground state in 
the S 3 A is a normal Fermi liquid. According to numeri- 
cal theories based on DMFT, 20 - 22 - 24 ) on the other hand, 
when U is sufficiently large the ground state appears to 
be a highly degenerate insulator for not only n = 1 but 
also n ~ 1. The assumption in the early Kondo lattice 
theory should be critically examined, because the ground 
state may be a highly degenerate insulator, as discussed 
above. However, whether or not the ground state in the 
S 3 A, DMFT, or DCPA can be a highly degenerate in- 
sulator itself is also an issue to be critically examined, 
because the third law is unlikely to be broken in a physi- 
cally relevant model or the Hubbard model with finite U. 
It is another issue whether or not the ground state can be 
a highly degenerate insulator beyond the S 3 A, DMFT, 
or DCPA. There is evidence that the ground state is a 
singlet 25 ) even within the constrained Hilbert space, 26 ) 
as discussed below. 

The Heisenberg model is derived from the half-filled 
Hubbard model. Fazekas and Anderson proposed that, in 
the two-dimensional triangular lattice, the ground state 
is a spin liquid stabilized by the formation of a local sin- 
glet or a resonating valence bond (RVB) on each pair of 
nearest neighbors by the superexchange interaction; 27 ) 
this spin liquid is called an RVB spin liquid and this 
stabilization mechanism is called an RVB stabilization 
mechanism or simply an RVB mechanism. The ground 
state of an RVB spin liquid is a singlet. 25 ) Although an 
RVB spin liquid is an insulator, it should be distinguished 
from an insulator whose ground state is highly degener- 
ate. In this paper, an insulator whose ground state is a 
singlet 25 ) is called a spin liquid and an insulator whose 
ground state is highly degenerate is only called a Mott 
insulator. 

The t-J model is also derived from the Hubbard model 
and is reduced to the Heisenberg model when n = 1. 
Anderson proposed that, in the two-dimensional square 
lattice, the ground state for n ~ I is an exotic electron 
liquid, 28 ) which is called an RVB electron liquid. The 
ground state within the constrained Hilbert space is cer- 



tain to be a singlet 25 ) for any n; it is stabilized by the 
RVB mechanism when n ~ 1 and is a normal Fermi liq- 
uid when n ~ 0. In general, no symmetry breaking occurs 
when a singlet ground state is transformed into another 
singlet ground state as a function of a parameter such as 
n, or the adiabatic continuation holds between the two 
sing lets. 29 ) Thus, it is an issue whether the ground state 
for n ~ 1 is really an exotic liquid or is simply a normal 
Fermi liquid. The RVB stabilization is due to the Fock- 
type term of the superexchange interaction, either when 
the ground state is an RVB electron liquid 30 ) or when 
it is a normal Fermi liquid. 31,32 ) When the RVB mech- 
anism is considered in the Hubbard model, the ground 
state is certain to be a singlet, 25 ) an exotic or normal 
Fermi liquid. 

The Bcthc ansatz solution for the Hubbard model in 
one dimension was obtained by Lieb and Wu. 33 ' Two 
chemical potentials arc defined by two processes of adding 
and removing a single electron: fi + (N) — E g (N + 1) — 
E g (N) and fi-(N) = E g (N) - E g (N - 1), where E g (N) 
is the ground-state energy for the electron number N. 
Then, a gap is defined by 

e 9 (AT) = MAO - MAO. (1-1) 

The number of unit cells is denoted by L. According to 
the Bethe ansatz solution, e g (L) > for any nonzero 
U in the thermodynamic limit of L — > +00. Then, Lieb 
and Wu argued that the half-filled ground state is an 
insulator for any nonzero U . If there is no other solution 
besides the Bethe ansatz solution, 34 ) the ground state for 
finite U is a singlet. 25 ) The insulator argued by Lieb and 
Wu, which is a singlet, 25 ) is a spin liquid but is not a 
Mott insulator, which is highly degenerate. 

The Bethe-ansatz solution is for the canonical ensem- 
ble, in which the electron number N is a constant of 
motion and is definitely an integer. In the grand canoni- 
cal ensemble or in the presence of an electron reservoir, 
on the other hand, N is not a constant of motion so 
that the quantum-mechanically averaged electron num- 
ber, which is denoted by (N), is not an integer, in gen- 
eral. Thus, whether the gap e g (N) is zero or nonzero 
is not a relevant criterion whether the ground state is 
a metal or an insulator. The compressibility is defined 
by k — (1/L)(d (N) /dfi), where /x is the conventional 
chemical potential. Since (N) is not an integer in gen- 
eral, what can be concluded from e g (L) > is that the 
half-filled ground state is either an insulator, in which k 
is exactly zero, or a metal in which k is infinitesimally 
small but nonzero in the thermodynamic limit. 35 ) Even 
if e g (L) > in the canonical ensemble, it is necessary to 
examine whether or not a gap opens in the single-particle 
excitation spectrum, which is defined in the grand canon- 
ical ensemble, of the half-filled ground state. 

The main purpose of this paper is to study the ground 
state of the Hubbard model in the presence of an elec- 
tron reservoir, but within the constrained Hilbert space 
where no order parameter exists. This paper is organized 
as follows: Preliminaries are given in § 2. The self-energy 
of electrons is decomposed into the single-site and mul- 
tisite self-energies. In § 3, it is proved that the single-site 
self-energy is that of a normal Fermi liquid even if the 
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multisite self-energy is anomalous, which leads to another 
proof that the ground state is a normal Fermi liquid in 
the S 3 A. An infinitesimally small reservoir effect is cru- 
cial in these proofs. In § 4, the nature of the ground state 
is studied and it is shown that the ground state cannot 
be a Mott insulator even beyond the S 3 A. In § 5, it is 
shown that the normal Fermi liquid in the S 3 A is fur- 
ther stabilized by the RVB mechanism. A discussion is 
given in § 6 and a conclusion is given in § 7. In Appendix 
A, an inequality is proved for the proofs in § 3. In Ap- 
pendix B, intersite exchange interactions are reviewed 
for the study in § 5. In Appendix C, it is argued based 
on the Curie- Weiss law for local moment magnetism that 
a Mott insulator is only possible as a high temperature 
phase. 

2. Preliminaries 

2.1 Electron reservoir 

The total Hamiltonian is composed of three terms: 

n = n a + n b + v. (2.1) 

The first term is the Hubbard model on a hypercubic 
lattice in d dimensions, which is called an a sublattice: 

f 



n " = e « E Ui ° ~ ~n E a h a 3°+ u Y n ^ ni ^ ( 2 - 2 ) 



(ij)a i 

where n ia — a\ a a iai e a is the band center, —t/Vd is the 
transfer integral between nearest neighbors (ij), and U 
is the on-site repulsion; here and in the following part 
of this paper, for convenience, the transfer integral is 
defined in such a way that it includes the dimensional 
factor d such as —t/y/d between nearest neighbors. 
The second term represents an electron reservoir on a 
hypercubic lattice in d dimensions, which is called a b 
sublattice: 

*6 



(2.3) 



where e& is the band center and —tb/Vd is the transfer 
integral between nearest neighbors. The third term rep- 
resents an infinitesimally small but random hybridization 
between the Hubbard model and the reservoir: 



V 



= * EE[< 



(2.4) 



where vu^ is a hybridization matrix between the zth site 
in the a sublattice and the jth site in the b sublattice in 
a set of 1Z: for convenience, an infinitesimally small but 
nonzero numerical constant 



A = ±0+, 



(2.5) 



is introduced. In order to recover the translational sym- 
metry in an ensemble averaged system, it is assumed that 

{ V (ij)} = and i V (ij) V (i'j')} = n hM 2 %)(i'j')' whcrC 

((•••)) represents an ensemble average over 1Z and is 
the density of pair of hybridization sites (ij) per unit 
cell referring to the a sublattice. The numbers of unit 
cells in the a and b sublattices are denoted by L and 
Lb- The thermodynamic limit is assumed: L — ► +00 and 
Lb/L — > +00. 



When the lattice constant of both the a and b sublat- 
tices is £, the structure factor of the sublattices is given 
by 

d 

/,(fc) = $^cos(M). ( 2 - 6 ) 

When [7 = 0, the Green function for electrons in the 
Hubbard model averaged over the ensemble is given by 

^ te -»> = le , + ,-^- 1 BW + iW (2 ' 7) 

where £; = (21 + ljirk^T, with I an integer, is the 
fermionic energy, \i is the chemical potential, 

E(k) = -2(t/Vd)f e (k), 



(2.8) 



and r(ie/) is the self-energy due to scatterings from the 
random hybridization or a reservoir term. The band- 
width of E(k) is W = 0(|i|). Since A = ±0+, the second- 
order perturbation is sufficiently accurate to treat the 
reservoir term, so that 



r(i £l ) = n h \ 2 \v\ 2 -^Yl 



E b (kY 



(2.9) 



where E b (k) = —2(tb/\/d)fe(k). It is assumed that 

Imr(+iO)<0. (2.10) 

The reservoir term ensures the quality of the grand 
canonical ensemble that the number of electrons in the 
Hubbard model, (N) = ^2 icr (ni a }, is a continuous func- 
tion of the chemical potential fx. It is assumed that in 
the half-filling case, in which (N) = L, there exists the 
particle-hole symmetry in the total system composed of 
the Hubbard model and the reservoir averaged over the 
ensemble. 

2.2 Fermi surface condition 

The Anderson model is defined by 



Ti-A = £d Y nda + E E c(k)c\ a c ha + Un^ridi 

(J ka 

+ J2 (V k c{ a d a + V^dUk*) , (2.11) 

where €d is the level of d electrons, n^a = d\d a , E c (k) is 
the dispersion relation of conduction electrons, Vk is the 
hybridization matrix, and La is the number of unit cells; 
for convenience, the on-site U is assumed to be the same 
as that in the Hubbard model. The Green function for d 
electrons is given by 

-1 



G a (m) 



isi+jl-ed- 



nj iei - e' 



(2.12) 

where fi is the chemical potential, E CT (ie;) is the self- 
energy, and A(e) is the hybridization energy defined by 

A(e) = -^^|V fe | 2 5[e + £-£ c (fc)]. (2.13) 
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The Fermi surface of conduction electrons, which is de- 
fined by E c (k) = fr, exists provided that 

A(0) > 0. (2.14) 

This condition is called the Fermi surface condition. 

The s-d model is another effective Hamiltonian to 
study the Kondo effect. According to Yosida's pertur- 
bation theory 36 ) and Wilson's renormalization group 
(RNG) theory, 37 ^ the ground state is a singlet or a nor- 
mal Fermi liquid except when J s -d = 0, where J s -d is the 
s-d exchange interaction. Since the s-d model is derived 
from the Anderson model, the result for the s-d model 
implies that the ground state of the Anderson model 
is also a normal Fermi liquid. According to the Bethe 
ansatz solution for the Anderson model with nonzero 
and constant A(e), the ground state is actually a normal 
Fermi liquid. 38 ~ 41 ) In general, the nature of the ground 
state depends only on relevant low-energy properties, 
such as A(0), and high-energy properties only quanti- 
tatively renormalize the ground state, as demonstrated 
by RNG theories for the s-d model. 37 ' 42 -* When the Fermi 
surface condition (2.14) is satisfied, the ground state of 
the Anderson model is a normal Fermi liquid except when 
U/A(0) = +oo and (rid) = 1, which corresponds to the 
case of J s -d = in the s-d model. 

3. Self-Consistent Single-Site Self-Energy 

In this section, no assumption is made for the on-site 
U, temperature T, and the electron density n — (N) /L. 
Within the constrained Hilbert space, in general, the 
Green function for electrons in the Hubbard model av- 
eraged over the ensemble is given in the wave-number 
representation by 

Ga{l£uk) = i £l + ^ - e a - E(k) - S ff (ie,, k) - r(ie,) ' 

(3.1) 

where S CT (ie/, fe) is the self-energy; it is given in the site 
representation by 



the single-site self-energy is given by the self-energy of 
the Anderson model: 



Rijofa) - }^e* (fl *-^ ) G CT (i £; ,fe), 



(3.2) 



where Ri is the position of the zth unit cell. 

Consider the Feynman diagrams for the self-energy of 
electrons in the site representation. Even if ensemble- 
averaged vertex corrections due to V appear in the di- 
agrams, they can be ignored because they are 0(A 4 ), 
where A = ±0 + . If only lines of the on-site U and the 
site-diagonal Green function, R iicr (iei), appear in a dia- 
gram, it is a single-site diagram. If a line of the site-off- 
diagonal Green function, Rij a {\e{) with i ^ j, appears, it 
is a multisite diagram. According to this classification the 
self-energy of electrons is decomposed into the single-site 
E CT (ie/) and the multisite AS CT (i£;,fc): 



£ CT (ie ( ,fc) = £ CT (ie ( ) + A£ CT (i£;,fc) 



(3.3) 



The same on-site U appears in the the Feynman dia- 
grams of the Hubbard and Anderson models. If — \i 
and A(e) are determined to satisfy 



E ff (i £i ) - E„(i£i). 
Equation (3.4) cannot be satisfied unless 

e a - A* = e d - fr- 
it follows from eq. (3.4) that 

A(e)=Im fe ff (e + i0)+i?^(e + i0) 



(3.5) 
(3.6a) 

(3.6b) 



Equation (3.4) or (3.6) is the condition that the calcula- 
tion of the single-site self-energy is mapped to a problem 
of self-consistently determining and solving the Ander- 
son model; in this paper, eq. (3.4) or (3.6) is called the 
mapping condition to the Anderson model. It should be 
noted that A(e) or the mapped Anderson model itself 
may depend on T. 

The mapping condition (3.6b) is iteratively treated to 
obtain the eventual self-consistent A(e); not only A(e) 
but also the single-site S CT (e) and the multisite AE^fe + 
iO, k) should be self-consistcntly determined with each 
other to satisfy the mapping condition (3.6b). It is proved 
in Appendix A that 



A(e) > -ImT(e-l-iO). 



(3.7) 



According to eqs. (2.10) and (3.7), the Fermi surface 
condition (2.14), A(0) > 0, is satisfied at each step of 
the iterative process, so that the ground state of the 
self-consistently determined or mapped Anderson model 
should be a normal Fermi liquid even if the multisite self- 
energy is anomalous or the ground state of the Hubbard 
model is a non-normal Fermi liquid. Thus, it is proved 
that the single-site self-energy for the ground state is that 
of the normal Fermi liquid in the Anderson model. Then, 
it is also proved that, when the multisite self-energy is ig- 
nored in the S 3 A, the ground state of the Hubbard model 
is a normal Fermi liquid. 

Consider the mapped Anderson model. Then, intro- 
duce an infinitesimally small Zeeman energy h — gfisH. 
Since the ground state of the Anderson model is a nor- 
mal Fermi liquid, the self-energy of the Anderson model 
is expanded in such a way that, for example, at T = K, 



E ff (e + iO) = So + (1 - c )e - (a + i~f)e 2 /k B T K 
-(l-4, s )^ah + 0[e 3 /(k B T K f] , 



(3.E 



(3.4) 



where So, 4> c > 0, <p s > 0, a, and 7 > are all real; 
7k or fceTk is the Kondo temperature, which is the en- 
ergy scale of local quantum spin fluctuations in not only 
the Anderson model but also the Hubbard model. Ac- 
cording to the proofs made above in this section, the 
expansion (3.8) is relevant even if the ground state of 
the Hubbard model is a non-normal Fermi liquid, and 
fesTk is definitely nonzero although it may be extremely 
or infinitesimally small. 

4. Nature of the Ground State 

In this section, the nature of the ground state in the 
weak and strong-coupling regimes, < U/\t\ < +00, 
is studied under only the assumption that the multisite 
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sclf-cnergy AS (J (e + i0,fc) is analytic in the upper half 
plane; it may be convergent or divergent on the real axis. 
The single-particle excitation spectrum is given by 

p(e) = --lmR lla (e + iO) = --ImG CT (e + iO), (4.1) 

IT IT 

where the mapping condition (3.4) is made use of. The 
density of states in the Hubbard model is the same as 
that in the mapped Anderson model. It follows from 
cqs. (3.8) and (4.1) that 



-— Im 

IT 



M{e)-E{k)- AE CT (e + iO,fc) 

(4.2a) 

n -1 



M(e) 



it J e - e' + iO 



(4.2b) 



where Eik) is defined by Eq. (2.8); when |e| <C k B T K , 

M(e) = Ap + 4> c e + {a + ij)e 2 /k B T K + O [e 3 /(k B T K ) 2 ] , 

(4.3) 



where 



Ap, = p - e a 



(4.4) 



The half-filling case (n = 1) is first studied. Since the 
particle-hole symmetry exists in not only the Hubbard 
model but also the Anderson model, it follows that 



%> = \v 

Ap, = 0, a 



p-e d = p 



0, A(e) = A(-e), 



and 



A(e') 



mtA(O). 



(4.5) 



(4.6) 



(4.7) 



e=0 



e - e' + iO 

Then, it follows from eq. (4.2b) that 

p(0) = 1/ [ttA(O)] . (4.8) 

When p(0) is vanishing or the ground state continuously 
approaches to an insulating state, A(0) is diverging such 
that A(0)/\t\ > and A(0)/\t\ -» +oo. Even if p(0)\t\ -» 
+ , the Fermi surface condition is definitely satisfied. 

In the S 3 A, AE CT (e + i0,fe) = 0. The ground state is 
definitely a normal Fermi liquid and no gap opens in p(e). 
Since Ap — 0, as shown in eq. (4.6), the Fermi surface is 
a hypersurface defined by 



E(k) = 0. 
It follows from eq. (4.2a) that 

p(°) = t£W0]- 



(4.9) 



(4.10) 



Thus, p(0) does not depend on U; A(0) does not depend 
on U either according to eqs. (4.8) and (4.10). 

Beyond the S 3 A, AE CT (e + iO, k) can be nonzero. As- 
sume that it is continuous and finite at e = for any k. 
According to the particle-hole symmetry, it follows that 

E{k) + RcA£ (T (+i0,fc) = 0, (4.11) 

for any k defined by eq. (4.9). Then, p(0) > 0, i.e., no 



gap can open in p(e). The ground state is a metal; it may 
be a normal or anomalous Fermi liquid. 

Since AE CT (e + iO, fe) is analytic, it can only diverge 
at a pole or an end-point of a cut; it cannot diverge on 
a line or an area. Then, assume that AS CT (e + iO, fe) is 
diverging as e — > ±0 on the real axis: 



lim |A£ CT (e + i0,fc)| = +oo. 



(4.12) 



If the divergence (4.12) occurs for a part of fe's defined by 
eq. (4.9), p(0) > and the ground state is a metal. If the 
divergence (4.12) occurs at least for all the fe's defined 
by eq. (4.9), p(0) = 0. Since i-fe 2 / (feB^k) is nonzero and 
AE (7 (e + i0, fe) is finite for e ^ 0, p(e) > for at least < 
|e| <C k B Tx- Then, a pseudogap opens and the ground 
state is a gapless semiconductor. Even if the divergence 
(4.12) occurs for any fe, the ground state cannot be an 
insulator with a complete gap. 

It is straightforward to extend the above analysis to 
nonhalf fillings (n/ 1). If no multisite self-energy is con- 
sidered in the S 3 A, the ground state is a normal Fermi 
liquid. If the multisite self-energy AE (J (e + iO, fe) is con- 
tinuous and finite at e = 0, the ground state is a metal. 
Even if the divergence (4.12) occurs for any fe, the ground 
state cannot be an insulator with a complete gap but can 
only be a gapless semiconductor. 

5. RVB Stabilization Mechanism 

In this section, the ground state in the strong-coupling 
regime, 1 <C U/\t\ < +oo, is studied. According to pre- 
vious studies for the t-J model, 31 ' 32 ' the Fermi liquid in 
S 3 A is further stabilized by the Fock-type term of the 
superexchange interaction. In this section, the study for 
the t-J model is extended to the Hubbard model. 

In the second-order perturbation in t, the superex- 
change interaction arises from the virtual exchange pro- 
cesses allowing empty and double occupancies. 43 ' In field 
theory, it arises from the virtual exchange of a pair exci- 
tation of electrons across the Hubbard gap: 44-46 ' 



J s (q) = 2(J/d)fe(q), 
where fe(q) is define by Eq. (2.6) and 
J Id = -4at 2 /dU, 



(5.1) 



(5.2) 



where a is a numerical constant; a — 1 when the band- 
widths of UHB and LHB are ignored, but < a < 1 
when they are considered. 47 ' An effective three-point 
single-site vertex function in spin channels is given by 
the expansion coefficient (f> s of the single-site self-energy, 
as shown in Appendix B. 

When the expansion (3.8) is used, the coherent part of 
the Green function is given by 

Gerfe, fe) = h • • • , 

(f> e e + Afi-E(k)- AE„(i£,,fc) 

(5.3) 

where Ap, is defined by eq. (4.4); (a + i-f)e 2 / (k B T K ) is 
ignored. When only the coherent part is considered, the 
Fock-type term is determined from 

A£ CT (k*, fe) = k B Tj ]T e^ 0+ j4> 2 J s (k - p) 
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i(f) c ei +A(jl- E(p) - AS CT (ie ; ,p) 
It follows from this self-consistent equation that 



AE ff (e + iO,fc) = - 



1 - 



cj(j/d)f t (k), 



(5.4) 



(5.5) 



where 



3 
d 



(0 s /0 o ) 2 i^0[-C(fc)/|t|]/,(fc), (5.6) 



where 0(e) is defined by 0(e > 0) = 1 and 0(e < 0) = 0; 
and £(fc) is the pole of eq. (5.3), so that 



£(fc) = ^ \-2(t*/Vd)Mk) - A J , (5.7) 



where 



t* =t- {l/8)4> e cj(J/Vd). (5.8) 

The Fock-type term is of higher order in 1 /d. 

Since the Fock-type term (5.5) does not depend on 
energy e so that it is normal, the ground state in this ap- 
proximation is a normal Fermi liquid. Then, £(fe) is sim- 
ply the dispersion relation of quasi-particle in the nor- 
mal Fermi liquid. According to the Fermi surface sum 
rule, 7,8 ) A/x can be determined by 



= 7 £0H;(*O/I*|]- 



According to the Fermi-liquid theory, 7 ' 8 -* the low- 
temperature specific heat is proportional to T such as 
C = IcT + • • • and the specific coefficient is given by 



7c 



Vfcl^O), 



(5.10) 



where 



P(0) 



i^<y[-2(t7v^)//(fc)-A/i]. (5.11) 



Since J/\t\ < 0, the Fock-type term enhances the band- 
width of quasi-particles, so that the Fermi-liquid ground 
state in the S 3 A is further stabilized by the RVB mecha- 
nism. However, p(0) is reduced by the RVB mechanism; 
A(0) is enhanced according to eq. (4.8) 

Low-energy local quantum spin fluctuations, whose en- 
ergy scale is fceTk, are renormalized by the superex- 
change interaction or high-energy intersite spin fluctu- 
ations, whose energy scale is 0{U). In such a situa- 
tion, the low-energy scale fceTk is simply proportional 
to the renormalized bandwidth of quasi-particles, i.e., 
fceTk oc \t*\. In particular, when <j> e 3> 1 or <f) c — > +oo, 



k B T K <x\J\. 



(5.12) 



When U is finite, &bTk i s definitely nonzero in this ap- 
proximation, even for the just half-filled case. 

Consider the half-filled case and assume that d is finite. 
In the limit of U/\t\ — > +oo with t 2 /U kept constant, the 
Hubbard model is reduced to the Heiscnberg model with 
the superexchange interaction constant 

J/d= -At 2 /{dU), (5.13) 

between nearest neighbors; this limit is called the Heiscn- 



berg limit. The probabilities of empty and double oc- 
cupancies in the Hubbard model are 0(t 2 /U 2 ), so that 
those in the mapped Anderson model are also 0(t 2 /U 2 ). 
In the Anderson model, the expansion coefficient <j> e of 
the self-energy is inversely proportional to the probabil- 
ities, so that <f> c — 0(U 2 /t 2 ). In the Heisenberg limit, it 
follows that t/t* — > or 

^^0, (5.14) 

(l/8)4> c cj(J/Vd) 

which means that the normal Fermi liquid is totally sta- 
bilized by the Fock-type term. It follows that 

p(0)|t| — 0, -ycocVd/\J\. (5.15) 

Note that the Fermi surface condition is definitely sat- 
isfied even in the limit of p(0)|t| — > 0, as shown in § 4. 
Since p(0) is vanishing but the T-linear specific-heat co- 
efficient 7c is nonzero and finite, a normal Fermi liquid 
in the Heisenberg limit is frustrated as much as an RVB 
spin liquid in the Heisenberg model. 27 ' It is plausible 
that the Fermi liquid in the Heisenberg limit is simply 
an RVB spin liquid in the Heisenberg model or the adi- 
abatic continuation 29 ) holds between the two liquids, at 
least in the approximation of this section, where only the 
Fock-type term is considered beyond the S 3 A. 



(5.9) 6. Discussion 



Suppose that the infinitesimally small reservoir term 
exactly vanishes such that T(e + iO) = 0. 48 ' In this case, 
the Fermi surface condition (2.14) may or may not be sat- 
isfied so that the ground state of the mapped Anderson 
model may or may not be a singlet. If the ground state 
of the Anderson model is a singlet, the ground state of 
the Hubbard model is a normal Fermi liquid in the S 3 A, 
as studied in § 3. Then, if the ground state of the Hub- 
bard model is not a normal Fermi liquid in the S 3 A, 
the ground state of the Anderson model is not a singlet. 
Thus, if a non-normal Fermi-liquid ground state is pos- 
sible in the S 3 A, it should be highly degenerate. This is 
the reason why only an insulator whose ground state is 
highly degenerate is called a Mott insulator in this paper. 

In the so called Hubbard I approximation, 2 ) for exam- 
ple, when n = 1 the self-energy is given by 



U U 2 



i0)' 



(6.1) 



which has a pole at e = — iO, so that the ground state is 
an insulator with a complete gap. In general, the ground 
state can be an insulator with a complete gap only if the 
self-energy £er(e + iO, ft) has a pole at e = — iO; the in- 
sulating ground state is highly degenerate, as discussed 
above, so that it is a Mott insulator. As proved in § 3, 
however, any self-energy with a pole at e = — iO cannot be 
a self-consistent S 3 A self-energy for satisfying the map- 
ping condition (3.6b) if once the reservoir term r(e + iO) 
is introduced. In the S 3 A, the ground state cannot be a 
Mott insulator if once the reservoir term is considered. 

Since the formulations of the S 3 A, DMFT, and DCPA 
themselves are exactly equivalent to each other, it is sur- 
prising that numerical theories based on DMFT 20 ' 21 ' 24 ) 
are inconsistent with the proof in §3. For example, Bulla 
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studied the half- filled (n = 1) ground state within the 
constrained Hilbcrt space by Wilson's numerical RNG 
method based on DMFT. 24 ) According to this numerical 
RNG theory, the ground state can be a Mott insulator 
for finite U, which is inconsistent with the proof in §3 
that ground state is a normal Fermi liquid for any finite 
U . A possible explanation for the inconsistency is that 
the reservoir term, which is infinitcsimally small, is con- 
sidered in the proof but no reservoir term is considered 
in the numerical RNG theory. Another possible explana- 
tion for the inconsistency is that numerical theories for 
T = K and T > K cannot treat a normal Fermi liq- 
uid with extremely or infinitesimally small fceTk because 
of limited numerical accuracy and nonzero temperature. 
When the proof in §3 is taken into account, the numer- 
ical RNG theory 24 ' for T = K simply means that, 
when U > W and n = 1, fee 7k is extremely or infinitesi- 
mally small in the S 3 A, DMFT, or DCPA. On the other 
hand, a Mott insulator is only possible as a high temper- 
ature phase, as discussed in Appendix C. When Tk is 
extremely or infinitesimally small, numerical theories for 
T > K cannot treat a crossover between a normal Fermi 
liquid at T < T K and a Mott insulator at T » T K . If an 
insulating phase appears at a low temperature T > K, 
the appearance means that Tk < T or A;bTk is extremely 
or infinitesimally small in the insulating phase. 

In this paper, an unrealistic electron reservoir is con- 
sidered in order to recover the translational symmetry in 
an ensemble averaged system, which simplifies the for- 
mulation as the periodic boundary condition, which is 
also unrealistic, simplifies the formulation. It is desirable 
to consider a more realistic electron reservoir. In general, 
a highly degenerate ground state is unstable against an 
infinitcsimally small perturbation, so that the third law 
of thermodynamics holds in a physically relevant model. 
Thus, it is unlikely that there is an S 3 A solution in which 
the third law is broken, even when a more realistic elec- 
tron reservoir is considered or no electron reservoir is 
considered. 

Even if a Mott insulator is possible in the S 3 A with 
T(e + iO) = assumed, it is unstable if once the 
RVB mechanism is considered beyond the S 3 A, as dis- 
cussed below. Single-site properties in the S 3 A such as 
A(e), which is the hybridization energy of the Ander- 
son model, So, <j) e , and <j> 3 , which are expansion coeffi- 
cients of the single-site self-energy, and so forth should 
be self-consistently calculated with the Fock-type term of 
the superexchange interaction. 49 ) First, we assume that, 
even when T(e + iO) = 0, <j) c is finite beyond the S 3 A, 
and then we examine the relevance of the assumption 
of finite <p e . The examination is exactly in parallel with 
the study in § 5. Provided that <j) c is finite, even if 4> c is 
extremely large, the Kondo temperature fceTk is renor- 
malized by the RVB mechanism, so that k B T^ > 0(\ J|) 
or fceTk = 0(|J|) even if r(e + iO) = is assumed. 
The nonzero /cbTk means that, when only the Fock-type 
term is considered beyond the S 3 A, a normal Fermi liq- 
uid, which is characterized by <j) c < +oo, is stable against 
a Mott insulator, which is characterized by <j) c = +oo or 
the self-energy with a pole at e = — iO. Thus, the assump- 
tion of finite 4> c must be relevant. Since k B T^ > 0(\ J\) 



or k B T K = 0(\J\) beyond the S 3 A, whether k B T K is 
nonzero, extremely small, infinitesimally small, or ex- 
actly zero within the S 3 A is never any crucial issue. Since 
the ground state cannot be a Mott insulator beyond the 
S 3 A, whether or not the ground state can be a Mott in- 
sulator within the S 3 A is never any crucial issue either. 

It is interesting to study if higher-order multisite terms 
than the Fock-type term can be so anomalous that the 
ground state is a gapless semiconductor. The divergence 
of the self-energy at the chemical potential implies that 
the ground state is highly or infinitely degenerate. Thus, 
the third law of thermodynamics appears to be incom- 
patible with the divergent self-energy. Therefore, the 
ground state within the constrained Hilbert space is un- 
likely to be a gapless semiconductor in one dimension 
and higher dimensions. 

It is straightforward to show that an anomalous term 
proportional to e In |e| exists in the multisite self-energy 
in one dimension; thus, the ground state in one dimen- 
sion is the so called Tomonaga-Luttinger liquid, which 
is not a normal Fermi liquid. Since the anomalous term 
is continuous and finite at e = 0, no gap opens. If the 
truth is that the self-energy is continuous and finite not 
only for n ^ 1 but also n = 1, the ground state within 
the constrained Hilbert space is a metal in one dimension 
and higher dimensions, even beyond the S 3 A. 

Suppose that an electron or hole is added onto the 
chemical potential in the half-filled Hubbard model. 
When N is not a constant of motion, the quantum- 
mechanically averaged number (N) is not an integer, in 
general. If a quasi-particle band exists near the chemical 
potential, a quantum process is possible in which only 
a small or infinitesimally small fraction of the added 
electron or hole enters the quasi-particle band and al- 
most the entire component enters the reservoir; thus, no 
gap opens in the single-particle excitation spectrum. It 
is trivial that if no quasi-particle band exists a gap opens 
in the spectrum. When TV is restricted within integers, 
on the other hand, the quantum process discussed above 
is impossible. The entire component of the added elec- 
tron or hole remains in the Hubbard model. It is trivial 
that, when U » W, the gap e g (L) defined by eq. (1.1) is 
nonzero such that e g (L) = U — 0(W) > in one dimen- 
sion and higher dimensions. When N — L and U 3> W, 
a gap as large as U inevitably opens in the spectrum of 
adding a single electron or hole, even if a quasi-particle 
band exists near the chemical potential. 

Observable physical properties or observables of an 
electron liquid are directly related with the pair exci- 
tation spectrum, which is bosonic. Although observables 
are not directly related with the single-particle excitation 
spectrum p(e), which is fermionic, certain observables 
can be described by p(e) according to the Fermi-liquid 
theory in the grand canonical ensemble, 7, 8 ) as shown in 
cq. (5.10). Since the virtual exchange processes of allow- 
ing empty and double occupancies are possible, a char- 
acteristic energy of low-energy pair excitations is 0(| J\) 
in the strong-coupling regime. Then, it is quite reason- 
able that the energy scale of 0{\ J\) appears in p(e). Since 
it is unlikely that observables in the canonical ensemble 
are different from those in the grand canonical ensemble, 
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the gap t g (L) in the canonical ensemble is not directly 
related with observables. It should be examined how the 
gap e g {L) in the canonical ensemble are related with ob- 
servables or what physical significance e g (L) has. 

The relevance of the assumption in the early Kondo 
lattice theory 16-18 ) has been confirmed in this paper; 
thus, the Kondo lattice theory is simply a perturbative 
theory starting from a normal Fermi liquid constructed 
in the S 3 A to study a normal or anomalous Fermi liq- 
uid and an ordered state in the whole Hilbert space. 
Not only the further stabilization of the normal Fermi 
liquid in the S 3 A by the RVB mechanism, 31 ' 32 ) which 
is also studied in § 5 of this paper, but also various 
anomalous Fermi-liquid properties can be explained by 
the Kondo lattice theory: paramagnon effect, 16 ) meta- 
magnetic transition or crossover, 50 ) and so forth. The 
Curie- Weiss (CW) law for itinerant magnetism is also 
an anomalous Fermi-liquid property; in a paramagnetic 
phase of an itinerant electron magnet, the static suscep- 
tibility x s (0, q) obeys the CW law for only q near an 
ordering wave number, but not for all q. In previous 
papers, 45 ' 51 ^ the CW law for itinerant magnetism has 
already been studied by the Kondo lattice theory. Many 
anomalous Fermi-liquid properties have been observed in 
the cuprate oxide superconductor at T > T c , where T c is 
superconducting critical temperature: resistivity almost 
linear in T, a pseudogap, and so forth. According to ref. 
52, resistivity can increase almost linearly in T because 
of critical antiferromagnctic fluctuations in two or highly 
anisotropic quasi-two dimensions. According to ref. 53, 
a pseudogap can open because of critical superconduct- 
ing fluctuations in two or highly anisotropic quasi-two 
dimensions. Such anomalous Fermi-liquid properties at 
T > T c can also be treated by the Kondo lattice the- 
ory when effects of antiferromagnctic or superconducting 
fluctuations are included in the multisite self-energy. 

An ordered state can also be studied by the Kondo lat- 
tice theory. In ref. 46, a theory of itinerant electron ferro- 
magnetism has already been developed based on a multi- 
band Hubbard model, in which the superexchange inter- 
action is ferromagnetic because of the Hund coupling. An 
early Fermi-liquid theory 54 ) of high-temperature super- 
conductivity, which was proposed just after the discov- 
ery of the cuprate oxide superconductor, 55 ) is consistent 
with the theory in this paper, so that it can be regarded 
as one of the simplest version of Kondo lattice theories 
of superconductivity. In a recent paper, 32 ) a theory of 
strong-coupling superconductivity has been developed to 
consider strong-coupling effects in the cuprate oxide su- 
perconductor. 

7. Conclusion 

In the canonical ensemble, the number of electrons is 
a constant of motion and is definitely an integer. In the 
grand canonical ensemble, on the other hand, the number 
of electrons within the Hubbard model is not a constant 
of motion because of the existence of an electron reser- 
voir, so that the quantum-mechanically averaged electron 
number is, in general, not an integer. Because of this dif- 
ference between the two ensembles, the single-particle 
excitation spectrum in the grand canonical ensemble, 



which is the spectrum of adding an infinitcsimally small 
number of electrons or holes, can be different from the 
spectrum of adding a single electron or hole in the canon- 
ical ensemble. According to the Fermi-liquid theory, the 
spectrum in the grand canonical ensemble is related to 
observable physical properties such as the specific-heat 
coefficient; therefore, the spectrum in the canonical en- 
semble is, in general, not directly related to observable 
physical properties, so that whether or not a gap opens 
in the spectrum in the canonical ensemble is not a rel- 
evant criterion whether the ground state is a metal or 
an insulator. Thus, the single-particle excitation spec- 
trum in the grand canonical ensemble has been studied 
to elucidate the nature of the ground state of the Hub- 
bard model or determine whether the ground state of the 
Hubbard model is a metal or an insulator, but within the 
constrained Hilbert space where no order parameter ex- 
ists. 

Within the constrained Hilbert space, the self-energy 
of electrons is decomposed into the single-site and mul- 
tisite self-energies, in general. It has been proved that 
the single-site self-energy for the ground state is that of 
a normal Fermi liquid even if the multisite self-energy 
is anomalous. If the multisite self-energy is continuous 
and finite at the chemical potential, no gap opens in 
the single-particle excitation spectrum; thus, the ground 
state is a metal, a normal or anomalous Fermi liquid. 
If the multisite self-energy is so anomalous that it is di- 
vergent at the chemical potential, only a pseudogap can 
open such that the ground state can be merely a gapless 
semiconductor; however, it is unlikely that the ground 
state is a gapless semiconductor since the third law of 
thermodynamics appears to contradict the divergence of 
the self-energy at the chemical potential. The ground 
state within the constrained Hilbert space cannot be an 
insulator with a complete gap nor a Mott insulator, at 
least when the on-site repulsion U is finite. When the 
on-site U is finite, a Mott insulator is only possible as a 
high temperature phase. 

In the supreme single-site approximation (S 3 A), in 
which the multisite self-energy is not considered, the 
ground state is a normal Fermi liquid. In the strong- 
coupling regime, the normal Fermi liquid is stabilized by 
the Kondo effect in the S 3 A and is further stabilized by 
the Fock-type term of the superexchange interaction be- 
yond the S 3 A. The stabilization mechanism is similar to 
that in the resonating-valence-bond (RVB) theory for the 
Heisenberg model. The normal Fermi liquid stabilized by 
the Fock-type term is frustrated as much as an RVB spin 
liquid in the Heisenberg model. It is plausible that, in two 
dimensions and higher, the adiabatic continuation holds 
between an RVB spin liquid in the Heisenberg model and 
a normal Fermi liquid in the Hubbard model. The Fermi 
liquid stabilized by the Fock-type term is a relevant un- 
perturbed state that can be used to study a normal or 
anomalous Fermi liquid and an ordered state in the whole 
Hilbert space by the Kondo lattice theory. 

Appendix A: Proof of the Inequality (3.7) 

An anomalous self-energy such as one given by 
eq. (6.1), which has a pole at e = — iO, can be used as a 
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trial input for the self-energy S CT (e + iO, fe) at the staring 
point of the iterative process according to the mapping 
condition (3.6b) to determine the Anderson model to be 
solved. Thus, it is only assumed in this Appendix that 
S CT (e + iO, k) is analytic in the upper half plane; it may 
have a pole or cut in the lower half plane and may be 
divergent on the real axis. Define the following real func- 
tions: 

5i(e,fe) =Rc[G; 1 (e + iO,fc)] , (A-l) 
S 2 (e,k) = lm[G- 1 (e + iO,k)], (A-2) 
S 2 (e) = -Im[r(e+iO) + t a (e+iO)] , (A-3) 

Yn{e) ^L^Sl(e,k) + sl(e,kY (A " 4) 



and 



and 



(A-5) 



L ^ Sl{e,k) + Sl(e,kY 
The site-diagonal Green function is given by 

Rua (e + iO) = Y 1 (e) - %Z X (e) . (A-6) 
According to the mapping condition (3.6b), 

3(6) 



A(e) = -Imr(e+iO) 



(A7) 



l?(e) + Z?(e)' 
where 

Z(e) = Z 1 (e)-S 2 (e)[Y 1 2 (e) + Z*(e)}. (A-8) 
In general, 

S 2 (e,k)>S 2 (e)>0, (A-9) 
for any k. It is trivial that Y (e) = Z (e), 

Y 2 (e) + Z 2 (e) = 1, (A- 10) 

and 

Zi(e) > S 2 (e)Y (e) = S 2 (e)Z (e). (A-ll) 
According to eqs. (A-10) and (A-ll), 

5(e) = ^(e) [Y 2 (e) + Z 2 (e)} S 2 (e)[Y?(e) + Z\(e)\ 
>S 2 (e) [-Y?(e) + Y (e)Y 2 (e)] 
+ S 2 (e) [-Zl(e) + Z (e)Z 2 (e)] . (A-12) 
Since inequalities of 



1^ [x + S^fc)] 2 



L^S?(e,ls)+S|(e,ls) 



and 



1 y> P + 



x + S 2 {e,k)Y 



) + Si(e,k) 



> 0, (A- 13) 



> 0, (A- 14) 



hold for any real x, i.e., Y (e)x 2 + 2Yi(e)x + Y 2 (e) > 
and Z Q (e)x 2 + 2Z 1 (e)x + Z 2 (e) > hold for any real x, 
it follows that 



-Z 2 (e) + Z (e)Z 2 (e) > 0. 



(A16) 



According to eqs. (A-9), (A-12), (A-15), and (A-16), it 
follows that 5(e) > 0. Thus, the inequality (3.7), A(e) > 
— Imr(e + iO), holds as a result of eq. (A-7) at each step 
of the iterative process. 

Appendix B: Intersite Exchange Interactions in 
Kondo Lattices 

It is assumed in this Appendix that U/W > 1. In gen- 
eral, the irreducible polarization function in spin chan- 
nels, which is denoted by 7r s (iw;, q), is also decomposed 
into the single-site 7r s (io;;) and the multisite A7r s (iw/, q): 

n s (iu u q) = n s (iwi) + Air s (iwi,q), (B-l) 

where u>i — Ilixk^I ', with I an integer, is the bosonic 
energy. The single-site 7r s (iw/) is also given by the polar- 
ization function of the self-consistently determined An- 
derson model. The spin susceptibilities of the Anderson 
and Hubbard models are given, respectively, by 

2jr s (iuJi) 



and 



Xs(iw;) 



Xs(iui,q) 



UTTsiiuJl)'' 



2ir s (\uj h q) 



(B-2) 



(B-3) 



. - [/7T s (iu> ( ,<?)' 

A physical picture for a Kondo lattice is that local spin 
fluctuations on different sites interact by an intersite ex- 
change interaction. According to this picture, the inter- 
site exchange interaction I s (iuoi , q) is defined by 



(B-4) 



1 - (l/4)J a (iwi,qf)x«(iwi)" 
It follows from eqs. (B-l)-(B-4) that 

I,(iui,q) = 2C/ 2 A7r s (i^,q){l + 0[l/L/x s (i^)]}- (B-5) 

When U /W > 1, terms of 0[1/Ux s (iu>i)] can be ignored. 

The exchange interaction I s (icoi , q) is composed of 
three terms: 45 - 46 - 51) 



I 



.(iL0i,q) = J s {q) + jQ{iuo h q)-AK{iuj h q). (B-6) 

The first term J s (q) is the superexchange interaction, 
which arises from the virtual exchange of a pair exci- 
tation of electrons across the Hubbard gap. The second 
term Jq (iw; , q) is an exchange interaction arising from 
the virtual exchange of a pair excitation of single-particle 
elementary excitations or quasi-particles near the chemi- 
cal potential in a normal or anomalous Fermi liquid. The 
third term is the remaining term, or the mode-mode cou- 
pling term among different modes of spin fluctuations, 
which corresponds to that in the self-consistent renor- 
malization (SCR) theory of spin fluctuations. 56 ' 

When the single-site irreducible three-point vertex 
function in spin channels is denoted by A s (k;,ie; + 
icji ;iu>i), it follows that 

A s (0,0;0) = ^[l-f/7r s (0)] 



-Y 2 (e)+Y Q (e)Y 2 (e) > 0, 



(A-15) 
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20s 



-{l + O[l/Ux.(0)]}, (B-7) 



UXs(O) 

according to the Ward relation. 57 ) When U/W > 
1, terms of O[l/Ux s (0)] can also be ignored. When 
eq. (B-7) is approximately used, the mutual interaction 
mediated by intersite spin fluctuations is given by 

i(C/A s ) 2 [ Xs (i^,q)-x s (i^)] = \4>lr s {\w u q), (B-8) 

where A s represents A s (0, 0; 0) and 



1 - (l/4)/ s (iw ( ,q)x s (iw;) 



(B-9) 



In eq. (B-8), the single-site term is subtracted. 

Multisitc or intersite terms can be perturbatively 
treated in terms of I s (\LOi,q) or I*{koi,q). The pcrtur- 
bative treatment is simply the Kondo lattice theory and 
is also the 1/d expansion theory. It should be noted that, 
according to eq. (B-8), <j) s is an effective three-point ver- 
tex function in this theory. 

Appendix C: Curie- Weiss law in a Mott Insula- 
tor 

When electrons are itinerant in a magnet, the Curie- 
Weiss (CW) law holds for only q near a magnetic or- 
dering wave number, which is the CW law for itinerant 
electron magnetism. When electrons are localized in a 
magnet, the CW law holds for any q, which is the CW 
law for local moment magnetism. The CW law for local 
moment magnetism can be used to define a Mott insula- 
tor, in which electrons are localized. 

First, suppose that d — ► +oo. Since magnetization is a 
leading-order effect in 1 /d, the CW law is also a leading- 
order effect in 1/d, in principle. In eq. (B-6), J s (q) and 
Jq{lo + iO, q) vanish for almost all q except for particu- 
lar q's such as high symmetric q's in the Brillouin zone 
and, if the ground state in the constrained Hilbert space 
is a normal Fermi liquid, nesting wave numbers in the 
Fermi surface; the particular q's are denoted by Q. It is 
obvious that — 4A(oj + iO, q) vanishes for all q. Thus, the 
susceptibility of the Hubbard model is such that 



Xs(uJ + iO,q) =Xs(u + iO), 
for almost all q except for Q and 

Xs (uJ + iO,Q) ^x s (uj + iO), 



(C-1) 



(C-2) 



for q = Q. The static susceptibility of the mapped An- 
derson model shows a crossover: 

Cl 



U0) 



k B T K ' 

C2 



T< T, 



K 



(C-3) 



T> T K 



fcB(T + c 3 T K )' 

where ci, Ci, and c 3 are numerical constants about unity; 
Tk may be extremely small. 

When kftT <C U, the T dependence of J S {Q) can be ig- 
nored. When T > T K , the T dependence of J Q (uj+iO, Q) 
can also be ignored. When Tk <C T <C U/ks, therefore, 
the static susceptibility Xs(0, q) of the Hubbard model 
obeys the CW law for all q, the CW law for local moment 
magnetism. If a Mott insulator is defined as an insulator 



that exhibits the CW law for local moment magnetism, it 
is only possible as a high temperature phase at T >• Tk 
or T > T K . 

When T < T K or T < T K , the T dependence of 
Xs(u + iO) is negligibly small; thus, that of Xs{u + iO, q) 
is also negligibly small except for Q. On the other hand, 
X«(w + iO,Q) can depend on T because Jq(0,Q) can 
depend on T. The precise T dependence of Jq(0,Q) or 
Xs (w + iO, Q) depends on the nature of the ground state. 
The static susceptibility Xs(0, <?) can obey the Curie- 
Weiss law, at least, if the ground state within the con- 
strained Hilbert space where no order parameter exists 
is a normal Fermi liquid. When the density of states for 
quasi-particlcs has a sharp peak near the chemical po- 
tential, the homogeneous one Xs(0, 0) obeys the CW law 
at T > Tc, where Tq is the Curie temperature. 51 ) When 
the Fermi surface of quasi-particles shows a sharp nest- 
ing at Qaf, the staggered one Xs(0) Qaf) obeys the CW 
law at T > 7n, where Tn is the Neel temperature. 45 ' The 
particular dependence of Xs(0, q) on q and T in infinite 
dimensions is a prototypic CW law for itinerant electron 
magnetism. It will be interesting to study how anoma- 
lous T dependence Xs( w + *0, Q) can have at T <C Tk 
or T < Tk, if the ground state within the constrained 
Hilbert space is not a normal Fermi liquid. 

In finite dimensions, according to the SCR theory of 
spin fluctuations, 56 - ) a CW like temperature dependence 
is possible due to the temperature dependence of the 
mode-mode coupling term — 4A(0, q). When p(e) has no 
sharp peak near the chemical potential and the Fermi 
surface shows no sharp nesting, — 4A(0, q) has the tem- 
perature dependence that gives the CW law for almost 
all q, but not for only q near a magnetic ordering wave 
number. This dependence on T and q is the so called CW 
law for local moment magnetism in an itinerant electron 
magnet. When pie) has a sharp peak near the chemi- 
cal potential or the Fermi surface shows a sharp nesting, 
— 4A(0, q) has a temperature dependence that suppresses 
the CW law. 51 ) This mechanism for the CW or anti-CW 
dependence is a higher order effect in 1/d. 
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